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Abstract 



The focus of this survey paper is on the distribution function 
FNp{t) for the largest eigenvalue in the finite N Gaussian Orthog- 
onal Ensemble (GOE, /3 = 1), the Gaussian Unitary Ensemble (GUE, 
O ' /5 = 2), and the Gaussian Symplectic Ensemble (GSE, /3 = 4) in the 

■ edge scaling limit of — > oo. These limiting distribution functions 

\ are expressible in terms of a particular Painleve II function. Compar- 

ts ' isons are made with finite simulations as well as a discussion of the 



universality of these distribution functions. 
1. Introduction 

In the well known Gaussian random matrix models of Wigner-Dyson- 



§ ! Mehta [O, O, O], the probability density that the eigenvalues lie in in 



finitesimal intervals about the points xi, . . . , xat is given by 

\ j<k 

where Cat/? is a normalization constant and 

1 for GOE, 
/?:=<( 2 for GUE, 
4 for GSE. 

We recall that for (3 = 1 the matrices are N x N real symmetric, for (3 = 2 
the matrices are N x N complex Hermitian, and for (3 = 4 the matrices 
are 2N x 2N self-dual Hermitian matrices. (For (3 = 4 each eigenvalue has 
multiplicity two.) 
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We are interested in 

En(3{0;J) := J ... JPNf3ixi,...,XN)dxi---dxN (1-1) 

Xj 

= probability no eigenvalues lie in J. 

For J = (t, oo), F]^fj{t) := i?jv/3(0; (t, oo)) is the distribution function for the 
largest eigenvalue, that is, 

Prob(Amax <t)= FNpit). 
If we introduce the n-point correlations 

N\ r r 

Rnpixi, . . . , Xn) '■= J-rz TT / ■ ■ ■ / Pn8{Xi, . . . , Xn) dXn+l " " " dXM , 

[N — n)\ Jn. Jn. 

and we denote by Xj the characteristic function of the set J, then we may 
rewrite ( p. . 1|) as 

EnI3{^]J) = / ■■■/ PNI3{xi,...,XN)T\{l-Xj{Xi))dXi---dXN 
JR JR ■ 

= ^2/3 X2) dxxdx2 - 

jjjjjj^3l3{Xl,X2,X3)dXidx2dX3 H (1.2) 

If we reinterpret P/v/? as the equilibrium Gibbs measure for like charges 
interacting with a logarithmic Coulomb potential (confined to the real line) 
subject to a harmonic confining potential, then everything we have said so 
far is valid for arbitrary inverse temperature (3 > 0. In this interpretation 
FN^it) is the probability, at inverse temperature P, that the interval (t,oo) 
is free from charge. 

2. Why (3 = 2 is the Simplest Case 



For P = 2 the n-point functions take a particularly simple form [|ri[ 

Rn2iXi, ...,Xn)= det {K^iXi, Xj) (2.1) 

where 

N-l 

KN{x,y) := Yl Vii^)Viiy) 

1=0 

(p{x)tlj{y) - ^(x)y(y) 
x-y 
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ipi{x) are the orthonormal harmonic oscillator wave functions, i.e. (pk{x) 

1 2 

Cfce~2^ Hk{x), Hk Hermite polynomials, and 



Using ( ^.11 ) in ( |1.2| ) we see that this expansion is the Fredholm expansion of 
the operator S whose kernel is KN{x,y)Xj{y); that is, 

EN2{0-J)=det{I-S). 

At this point we can use the general theory |jl6|, [l^ for Fredholm deter- 
minants of operators K whose kernel is of the form 



^{x)ij{y) -ij{x)'^{y) 
x-y 

where the ip and ip are assumed to satisfy 



dx \ip^ 



n{x) 



with Q{x) a 2 X 2 matrix, trace zero, with rational entries in x. It is shown 
in |13 that if 



then 



_d_ 

da. 



J := [j ia2j-i,a2j) 



logdet (I-K), j = l,...,2m, 



are expressible polynomially in terms of solutions to a total system of partial 
differential equations (a^ are the independent variables). 

In the finite GUE case 



n{x) 



-X V2N 
^ X 



and the theory gives for J = {t, oo) 

FN2i't) = exp {~J R{x) dx 

where R satisfies {' = d/dt) 

{R"f + A{R'f{R! + 2N) - A{tR' - Rf = 0. 
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This last differential equation is the a version of Painleve IV {Piv) 0; P^] - 
What is remarkable about this result is that the size of the matrix, N , (equiv- 
alently the dimension of the integral ( |1 . 1|) ) enters only as a coefficient in the 
above second order equation. 

3. Edge Scaling Limit 

The famous Wigner Semicircle Law states that if Pn{x) is the density of 
eigenvalues in any of the three Gaussian ensembles, then 

1- 1 /r, fTT \ \ - Vl — \x\ < 1, 

hm =pN[2aJNpx) = { L / 

/ Ar_ V V / .T > 1 . 



^p^'^2(T^Np V V ^ y |a;| > 1. 

Here a, a/v^, cr/^/T. (for [3 = 1,2,4, respectively) is the standard deviation 
of the Gaussian distribution in the off-diagonal elements and 

r N,f3 = i, 

Nf,= \ N,(3 = 2, 

{ 2N + 1, (3 = A. 

For the normalization here and in [^, 

a = I/V2. 

Perhaps less well known is that the following is also true 



FNp[2aJNp + x 



if a; < 

1 if a; > 



as — s> 00. The edge scaling variable, s, gives the scale on which to study 
fluctuations 0, |^, [l^ 

t = 2av^+-^. (3.1) 
The edge scaling limit is the limit N 00, s fixed, and in this limit |TB|, 



FN2{t) exp ^ (2^ - ■5)q'(a;)^ dx^ =: ^2(5) 
where q is the solution to the Pu equation 

q" = sq + 2q^ (3.2) 

satisfying the condition 

q{s) ~ Ai(s) as s — 00, (3.3) 
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with Ai the Airy function. 

We note that F2{s) is the P// r-function p, 0. There is a rather complete 
description of the one-parameter family of solutions to Pn satisfying the 
condition 

q{s; A) ~ AAi(s) 

as s — > oo and an analysis of the corresponding connection problem for 
s — > — oo (see 0, |^ and references therein). 

4. Cases P = 1 and 4 

In Dyson showed for the circular ensembles with /5 = 1 or 4 that the 
n-point correlations could be written as 

{Rn/3{Xl, . . . , Xn)f = det {^KN/jiXi, Xj)\lj^-^ 

where now 

Kfyfj3{x, y) = 2 X 2 matrix. 



Mehta generalized this result to the finite Gaussian ensembles and 
Mahoux and Mehta [10| gave a general method for invariant matrix models 
for both [3 = 1 and (3 = 4. Mehta's result implies that for /5 = 1 or 4 that 

{Er,p{0;J)f = det {I-Knp) 

where A^^r^ is an operator with 2x2 matrix kernel, or equivalently a 2 x 2 
matrix with operator entries. Explicitly for (3 = 1 (and A^ even): 



A'tVI = X ; 



/ S + i)®eLp SD-ilj®ip 

\ eS — e + eip ® Eif S + ecp ^ip 
where S, ip and ip are as before for /3 = 2, and 



Xj, 



ef{x) := / e{x -y)f{y)dy, 
Dfix) := 



dx 

with 



I if X > 0, 



e{x) := < ifx = 0, 
[-| ifx<0. 

We have used the notation a 6 for the operator with kernel a(x) b{y). 
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For (3 = A, (^4(0; J/v^)) is again a Fredholm determinant with 

/ S + ilj®eip SD -tfj^ip 
\ eS + eip ^ Eip S + Eip^ip 



K 



N4 
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In these Fredholm determinants are related to integrable systems for gen- 
eral J. Although |jl8[ treats exclusively the Gaussian ensembles, the methods 
appear quite general and should apply to other ensembles as well. 

5. Idea of Proof and Results 



Our derivation [|I^ rests on the fact that operator determinants may be 
manipulated much as scalar determinants, as long as one exercises some care. 
We first write 



K 



I XjD 



Nl 



{Se -Sip® e^lj)Xj {S + e(p(g) ip)Xj 

V Xj J \ {Se-e-eip®etlj)Xj {S + eip®-^)Xj 
This uses the trivial fact that De = I and the commutator identities 

[S, D] = ip ® ifj + ip ®(p, [e, S] = -Eif ® sip - eip ® e(p. 

By the general identity det(/ — AB) = det(/ — BA) the determinant is 
unchanged if the factors are interchanged and so we may work instead with 

{Se-eip®etlj)XjD {S + sip ® iIj)Xj 

{Se - e - e^®ei))XjD {S + £^®i))Xj 

Applying a pair of row and column operations, which is justified by the 
identity det(/ — A) = det(/ — BAB~^), we reduce this to 

{Se -e(f® e^)XjD + {S + e(f ® ^)Xj {S + eip ® ^)Xj 

-eXjD 

More row and column operations, which this time use det(/ — A){I — B) = 
det(J — A) det(/ — B), followed by factoring out 

I-SXj 

show that Eni{0; jy equals En2{0] J) times a determinant of the general 
form 

n 

det (/ - ^ Ofc ® Pk), 

k=l 
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whose value equals that of the scalar determinant 

det [6j^k-{aj,f3k))^^^^^^ ^^. 

Thus the evaluation of Eni{0; J) is reduced to the evaluation of certain inner 
products and the result for _E'jV2(0; J). 

To evaluate the inner products differential equations are derived for them. 
For J = [t, oo) and in the edge scaling limit these differential equations can 
be solved with the result that 

/'OO , , 

^xv-y ^zv^/- = '^"^ (5.1) 

where we recall 

-^2(5) = exp ^— y {x — s) q{xY dx^ . (5.2) 

Similary for /3 = 4 

F,{s/V2f = F2{s) i i . (5.3) 

We note that there are no adjustable parameters in Ff^. 

Using the known asymptotics of q{s) as s — ^> ±00, it is straighforward to 
solve (|3.2| ) numerically and produce accurate numerical values for Fp^s) and 
the corresponding densities fpis) = dFp/ds. The densities for the largest 
eigenvalue in each of the three ensembles in the edge scaling limit are dis- 
played in Fig. 1. Observe that for higher "temperature" (lower (3) the vari- 
ance of //^(s) increases as one would expect from the Coulomb gas interpre- 
tation. In Table 1 we give the mean (/i/3), standard deviation (cr/3), skewness 
{Sp), and the kurtosis {Kp) of the densities //3.Q 
6. How "Useful" are these Scaling Functions? 

Using a random number generator to simulate a real symmetric matrix 
of size X in the GOE, we calculate its largest eigenvalue Amax- For 
large the expected value should be approximately 



E (A^ax) = 2aViV + 



Ari/6 



where /xi is the mean of /i (see Table 1). In one such simulation (with a = 
I/V2) with = 500 and 5000 trials the mean largest eigenvalue is 31.3062 
which should be compared with the theoretical prediction of 31.353. 

""^The skewness of a density / is / (^^)'^ f{x) dx and the kurtosis is J (^^)^ f{x) dx—3 
where the —3 term makes the value zero for the normal distribution. 
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Probability densities 
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Figure 1: The probability densities //3(s) = dFp/ds, (3 = 1,2,4, for the 
position of the largest eigenvalue in the edge scaling limit. 



p 


/i/3 3/3 K/3 


1 


-1.20653 1.2680 0.293 0.165 


2 


-1.77109 0.9018 0.224 0.093 


4 


-2.30688 0.7195 0.166 0.050 



Table 1: The mean (/i/j), the standard deviation (ci/j), the skewness (Sjs) and 
the kurtosis (Kp) for the the densities 



For each largest eigenvalue Amax we compute a scaled largest eigenvalue 
s from 

/ — as 
^-ax = 2aVN + ^ 

and form a histogram of s values. In Fig. 2 we compare this histogram with 
the limiting density fi{s). In Table II we give the mean, standard deviation, 
skewness and kurtosis of this same scaled largest eigenvalue data. This should 
be compared with the limiting theoretical values in the first row of Table I. 
Similarly in Fig. 3 we give comparisons of finite A^(= 100) simulations in GUE 
and GSE with the limiting densities /2 and /4, respectively. The descriptive 
statistics of this data are given in Table II. 

7. Questions of Universality 

Recall that the finite N GOE is the unique^ measure on the space of 

^Uniqueness is up to centering and choice of standard deviation on, say, the ofF-diagonal 
matrix elements. 
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GOE Largest Eigenvalue Simulations 
.41 




s 



Figure 2: Histogram of scaled largest eigenvalues for N — 500, GOE for 5000 
trials. Solid curve is the limiting density fi{s). 













1 


-1.2(311 


1.2182 


0.273 


0.171 


2 


-1.7772 


0.8952 


0.133 


-0.089 


4 


-2.3154 


0.7057 


0.171 


0.023 



Table 2: The mean (///?), the standard deviation {ap), the skewness {S13) and 
the kurtosis {Kjs) for the scaled largest eigenvalue data from simulations in 
GOE {N = 500), GUE {N = 100) and GSE {N = 100). In each ensemble 
there were 5000 trials. 
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Largest Eigenvalue Simulations 




Largest Eigenvalue Simulations 




Figure 3: Histograms for the scaled largest eigenvalues in GUE and GSE 
each with N — 100 and 5000 trials. Sohd curves are the limiting densities 
/2(s) and /4(s), respectively. 
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N X N real symmetric matrices A that is both invariant under all orthogonal 
transformations and the matrix elements Aj^ with j > k are statistically 
independent. Similar statements hold for both the GUE and GSE. 

It is of interest to inquire to what extent the various scaling limits of 
these finite Gaussian ensembles are universal. For the edge scaling limit, 
the problem is to characterize the domain of attraction of the limiting laws 
(3 = 1,2,4, as given by (|5.1|) - (|5.3|) . The problem as so stated is open 
but there has been recent progress in showing that this domain of attraction 
contains more than the Gaussian case. The most studied case is that of 
invariant random matrix models of Hermitian matrices where the measure is 
of the form 

I^NidA) = Z^^ exp {-Tr{V{A)) dA 
and is a polynomial. The Gaussian case corresponds, of course, to V{x) = 

In the cases V{x) = and V{x) = j^x^ are analyzed. The nonuni- 
veral part of the problem is the transformation to the scaling variable s. For 
example, in the quartic case ( p.l|) is replaced by 

But once this is done, it follows from the results in |^ that the distribution 
of the largest eigenvalue converges in the edge scaling limit, as defined by 
the above equation, to the Gaussian result F2{s). 

In PI there is a an account of the Hermitian matrix model with measure 
of the form 

IJN{dA) = Z]^^exp{-NTT{V{A)) dA, 

with 

V{x) = ^x^ + ^x\ t<0,g>0. 

Note the insertion of the factor into the exponential. This has the ef- 
fect that the density of eigenvalues Pn{z) converges without any rescaling of 
variables to a limiting density p{z) given by 



27r 

with 

^1,2 
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In this set up the edge of the spectrum remains bounded and the transfor- 
mation to the edge scahng variable s is now 



It is proved that the distribution of the largest eigenvalue converges to 
the Gaussian result F2{s). 

Physicists (see and references therein) have heuristic arguments that 
suggest that the domain of attraction of the invariant matrix models for 
Hermitian matrices contains all those potentials V for which the density 
vanishes at the edge of the spectrum as does the Wigner semi-circle, i.e. a 
square root. By tuning the potential one obtains densities that vanish faster 
than a square root — these models will be in a different universality class with 
regards to the edge scaling limit. 
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